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ABSTRACT 

The relation of the solar surface magnetic field with mesogranular cells is 
studied using high spatial 100 km) and temporal (~ 30 sec) resolution data 
obtained with the IMaX instrument aboard SUNRISE. First, mesogranular cells 
are identified using Lagrange tracers (corks) based on horizontal velocity fields 
obtained through Local Correlation Tracking. After ^ 20 min of integration, the 
tracers delineate a sharp mesogranular network with lanes of width below about 
280 km. The preferential location of magnetic elements in mesogranular cells is 
tested quantitatively. Roughly 85% of pixels with magnetic field higher than 100 
G are located in the near neighborhood of mesogranular lanes. Magnetic fiux is 
therefore concentrated in mesogranular lanes rather than intergranular ones. 

Secondly, magnetic field extrapolations are performed to obtain field lines 
anchored in the observed flux elements. This analysis, therefore, is independent 
of the horizontal flows determined in the flrst part. A probability density function 
(PDF) is calculated for the distribution of distances between the footpoints of 
individual magnetic held lines. The PDF has an exponential shape at scales 
between 1 and 10 Mm, with a constant characteristic decay distance, indicating 
the absence of preferred convection scales in the mesogranular range. Our results 
support the view that mesogranulation is not an intrinsic convective scale (in 
the sense that it is not a primary energy-injection scale of solar convection), but 
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also give quantitative confirmation that, nevertheless, the magnetic elements are 
preferentially found along mesogranular lanes. 

Subject headings: quiet Sun — Sun: activity — solar convection. 



Introduction 



Mesogranulation was historically introduced as a prominent scale imprinted on the hor- 
izontal photospheric flows calculated through local correlation tracking (LCT) of intensity 



images ( 
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Roudier et al.||1998, 


1999 



1988 1991 MuUer et al. 



positive and negative divergence of that flow as well as the time evolution of Lagrange tracers 
moving in it revealed cells with sizes between, say, 5 and 10 arcsec. Much debate ensued 
concerning whether the mesogranular flow patterns correspond to actual convection cells in 
that range of sizes rather than, e.g., to simple granule associations which persist in time 
(see, e.g. ICattaneo et al.||2001t iRoudier et al. 1120031 iRoudier & Muller||2004l iNordlund et al. 



2009 Matloch et al. 2009, 2010). Independent hints for the existence of a convective flow 
operating on those scales are therefore of importance. The study of the surface distribution 
of magnetic elements can provide such hints; as a minimum, it can constitute an alternative 
avenue, independent of the inaccuracies of LCT methods, to determine the properties of 



the mesogranular patterns. Such an approach has been employed by Dommguez Cerdeiia 



et al. (2003); Dommguez Cerdefia (2003); Sanchez Almeida (2003) using ground-based data 



and by Roudier et al. (2009) and Ishikawa & Tsuneta (2010) using satellite observations. In 



those papers, visual evidence was obtained that there is an association between magnetic flux 
structures (flux elements, transient horizontal fields) and the mesogranular pattern obtained 
through the study of horizontal flows. Yet, detailed quantitative studies and statistics that 
could put such an association on a firmer basis are still missing. 

The aim of this letter is to obtain quantitative information of the relation between 
photospheric magnetic flux distributions and mesogranular scales by using the observations of 



unprecedented quality provided by the Imaging Magnetograph eXperiment (IMaX; Martinez 



Fillet et al.||2010[ ) aboard the SUNRISE balloon-borne observatory ( [Barthol et al.||2010 
Solanki et al.|20l"0 ). IMaX provides time series of virtually seeing-free high spatial resolution 
(0.15 arcsec) images and magnetograms that constitute an ideal data set for the study of 
photospheric magnetism. Using the Sunrise/IMaX data, we combine information from (a) 
the velocity field gained through Local Correlation Tracking of intensity images; (b) the 
spatial patterns provided by the magnetograms; and (c) the field line structure obtained 
through extrapolation of the magnetogram data, to gain quantitative information concerning 
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patterns at intermediate scales between granulation and supergranulation. 



2. Data 

For this study we use sequences of images recorded with IMaX near the solar disk center 
on 2009 June 9. Images were taken at five wavelengths along the profile of the magnetic- 
sensitive Fel 5250.2 A line located at ±80 mA, ±40 mA from line center, and continuum at 
+227 mA. The estimated circular polarization noise is 5 x 10~^ in units of the continuum 
wavelength for non-reconstructed data and 3 times larger for the reconstructed one. IMaX 



has a spectral resolution of 85 mA and a spatial resolution of ~ 100 km (Martinez Fillet 
et al.||20l"0 ). The reduction procedure produces time series of images with a cadence of 33.25 
s, spatial sampling of 39.9 km and a field-of-view (FOV) of 32 x 32 Mm^. We use two time 
series, the first one comprising 42 snapshots (23 minutes) and the second one 58 snapshots 
(32 minutes). Magnetograms are derived from inversions of the observed Stokes parameters 



using the SIR code (Ruiz Cobo & del Toro Iniesta 1992). We call these magnetograms 



the reconstructed data. In this letter we mainly use the reconstructed data except in the 
last section (Sec. [s]) where, in addition, we use so-called calibrated data, i.e., magnetograms 
obtained using a proportionality law between Stokes-V (from non-reconstructed data) and 
Bz (for details of this method, see Martinez Fillet et al.||2010 ). The intensity maps are taken 



from the reconstructed data. 

Figure [l] displays a map of normalized continuum intensity (top) and the corresponding 
longitudinal magnetogram (bottom), the latter showing many internetwork fiux concentra- 
tions alongside stronger fiux elements probably belonging to the network. 



3. Mesogranulation and horizontal flows 

Following the traditional method, mesogranular patterns are outlined by using Lagrange 
tracers (popularly known as corks) arranged at time t = uniformly in a 2D grid that co- 
incides with the observational pixel matrix. Then, the corks are advected following the 
instantaneous horizontal velocity field along the whole duration of the time series (here we 
use the second time series: 32 minutes). The horizontal velocity field is determined using 



the local correlation tracking (LCT) algorithm by Welsh et al. (2004) applied to consecutive 
continuum intensity images. The correlation is performed in local windows weighted by a 
Gaussian function with FWHM=320 km. Often, in the literature, the velocity field deter- 
mined through LCT and used to advect corks is a time average of the noisy determinations 
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Fig. 1. — Top: normalized continuum intensity near the spectral line Fe I 5250.2 A. Bottom: 
Vertical component of the magnetic field vector [Bz) retrieved from inversions and clipped 
at ±50 G. The maps are taken from the first observed time series. 
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of the instantaneous velocity fields for consecutive snapshots (for a discussion about the 



effect of time averaging see e.g. Rieutord et al. 2000, 2001). The high signal-to-noise ratio 



the high spatial and temporal resolution and the absence of atmospheric distortion in the 
IMaX data make it possible to evaluate clean velocity fields and consequently to compute 
cork advection without time averaging. 



In Fig. |2]we show (white dots) the distribution of corks at t = 28 min, i.e., toward the 
end of the series, on the background of the continuum intensity image at that time. The 
corks clearly delineate cells of mesogranular size containing several granules. In fact, the 
vast majority of corks are located in the lanes between those cells (which we will call here 
mesogranular lanes or mesolanes for short). A clear mesogranular pattern with the majority 
of the corks concentrated in mesolanes is already visible after only some 20 minutes from 



t = 0. This is less than the time (2-4 hours) reported by Roudier et al. (2009) using the 
Hinode/SOT NarrowBand Filter Imager. 

In order to obtain quantitative measures for the mesogranular network we define a cork 
density function, Pcork, by counting the number of corks in each pixel in the image. To 
provide an upper bound for the width of the mesolanes, we scan the image vertically using 
whole horizontal lines and calculate the width (1/e of the maximum) of the peaks of the 
Pcork function in each line. We then derive an average peak width for the whole scan. If we 
use all peaks above pcork = 10, then the resulting upper bound for the lane width is 285 km. 
The largest concentrations of corks are located in narrow mesolanes: using all peaks above 
Pcork = 30, for instance, yields a width of 180 km (4.5 pixels). Given the random orientation 
of the mesolanes, their actual width is certainly below those numbers. Using a horizontal 
scan instead of a vertical one yields basically the same numbers (the maximum deviation is 
7%). 



4. Correlation between mesogranular lanes and field concentrations 



A first test of the relation between mesogranules and the surface magnetic field can 
be obtained by studying the spatial association between mesogranular cells and magnetic 
flux concentrations. Through visual inspection, various researchers have obtained indica- 
tions that the magnetic field concentrations with higher flux density are distributed at the 



boundaries of cells of mesogranular size (5 to 10 arc sec; see Dommguez Cerdena et al. 2003 



and Sanchez Almeida 2003); in fact, Dommguez Cerdena (2003) shows the preference of 



elements with flux density above 60 G to be located in regions of high negative divergence 



of the horizontal vlocity. Further visual evidence is provided by Roudier et al. (2009) who 



plot Stokes- V images obtained with Hinode/NFI on top of cork distributions showing a rela- 
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tionship between magnetic field concentrations and cork lanes (see also de Wijn et al. 2005 



Solanki et al. 2010) 



The conclusions at that level can be reinforced through our high-resolution maps by 
drawing contours of \Bz\ down to small values, like 30 G. The red and yellow contours in 
Fig. [2] correspond to flux density values of 30 and 50 Gauss, respectively. The figure exhibits 
flux concentrations mostly located at mesolanes. One can notice that flux concentrations 
with higher flux density (yellow contours: 50 G) correlate better with mesolanes. 

Going one step further, we can obtain quantitative estimates to test the visual impression 
gained by combining cork images and magnetograms. We use as a proxy for the mesolanes 
the locations with cork density above a threshold, > p^g, and use the values p^g = 2, 4, 8, 16 
corks / pixel. All of those values yield clear mesogranular lanes, as apparent, e.g., in the 
top-right and bottom-left panel of Fig. |3| drawn for p^g = 8 and 2, respectively. Increasing 
PcQ decreases the level of connectivity of the lanes. We then focus on the mesolanes defined 
by Pc > Pco ~ ^ study which fraction of the magnetic elements are really in (or near) 
them. To obtain a quantitative estimate, for each fixed threshold intensity l-Bzl^/j in the range 
(0, 400) G, we consider the set of pixels with B > l-B^I^/j and calculate which fraction of the 
set is located in 3-pixel neighborhoods of the mesolanes. The result is shown in the top left 
panel of Fig. [3} red curve. The curve reaches an approximate horizontal asymptote at the 
85% level for l-Bzl^/j above 100 G. The other curves in the figure contain the same results 
but using p^g = 2 (black curve), 4 (green) and 16 (orange). We see approximate horizontal 
asymptotes in all cases, and all are reached for \Bz\ff^ between 80 to 100 G. We conclude 
that the vast majority of pixel elements with B > 100 G is located in the near neighborhood 
of locations with high cork density. For the least restrictive case (p^g = 2), virtually all 
magnetic elements with B > lOOG are located in the neighborhood of the mesolanes. For 
large p^g (e.g., p^g = 16 c/p, orange curve), the resulting areas have less connectivity and 
do not delineate so clearly the mesogranular network. Correspondingly, the asymptotic 
percentage values for B > 100 G become smaller. The width of the 3-pixel neighborhoods 
of the mesolanes is shown in the bottom-right panel for the p^g = 8 corks/pixel case. The 
panel shows a fully developed network of mesolanes covering only 17% of the whole surface 
(which, of course, coincides with the l-B^I^^ — limit of the red curve). As an aside note, 
all curves remain basically the same if we remove the prominent network element located in 
the lower part of Fig. [2j 
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Fig. 2. — Cork distribution (white dots) outlining mesogranular cells at t=28 min. The 
background image is continuum intensity. Color lines: \Bz\ contours for 30 G (red) and 50 
G (yellow). 
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Fig. 3. — Top-left: percentage of pixels with \Bz\th above the value given in abscissas located 
within 3 pixels (~120 Km) from mesolanes. Mesolanes are defined through > 2,4,8,16 
corks/pixel for the black, green, red and orange curves, respectively. Top-right and bottom- 
left panels: distribution of pixels with > 8 (top right) and > 2 (bottom left) corks/pixel, 
respectively. Bottom-right panel: 3-pixel neighborhood of the points shown in the top-right 
panel. The calculations are done for similar time as Fig. [2] (t=28 min). 
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5. The separation between the footpoints of extrapolated field lines 

5.1. Method 



As a further test of the relation between surface fields and convection at different scales, 
in this section we study the statistics of footpoint separation between field lines linking the 
magnetic elements observed with IMaX. To that end, we calculate the magnetic field vector 
in the atmosphere using force-free field extrapolations (FFF) from the IMaX data. We use 
a code with weighted optimization method (Wiegelmann 2004 see also Seehafer||1978 ). The 
force-free assumption is equivalent to assuming that the electrical current j and the magnetic 
field B are parallel, or, using Ampere's law: 



V X B 



4tt . 



aB ; 



a in Eq. [T] measures the level of field line twist. The field values obtained with the IMaX data 
correspond to the height where the Fe I 5250.2 A line is formed. Although the photosphere 
is not a small-/3 plasma (with /3 the ratio of gas to magnetic pressures), there are theoretical 
and observational indications (Wiegelmann et al. 2010a||b Martinez Gonzalez et al. 2010) 
that the force-free assumption is acceptable when calculating extrapolations starting from 
those heights. We will be using a = 0.1/ L, with L x L being the IMaX field of view. In any 
case, Wiegelmann et al. (2010b), using the same IMaX data, show that values of a in the 



range {—4:/L,4:/L) lead to similar values in the statistical properties they analyzed for the 
extrapolated field lines. 

For computational reasons, we have calculated the extrapolations in a box of 389 x 389 
X 389 equally spaced grid points, keeping the original horizontal dimensions, leading to a 
cell size of ~ 80 km. Yet, test calculations with the original number of grid points (778 in 
each direction) yield essentially the same results. 

An example of magnetic field lines is shown in Fig. |4} A map of the vertical magnetic 
field component is shown at the bottom of the box corresponding to a region of about 10 x 10 
Mm^ centered on the point {x,y) = (10, 20) Mm in Fig. [!} The field line footpoints are seen 
to cluster on field concentrations in the domain. In agreement with the conclusions of 



Wiegelmann et al. (2010b), this figure shows a rough proportionality between the footpoint 



distance and the height reached by the field lines. 
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5.2. Statistics of footpoint separation 

Calling X the distance between footpoints, we calculate the probability density function 
(PDF) for X using as statistical ensemble for each time series the field lines in each snapshot 
and the whole collection of snapshots in the series. To include a field line in the study, 
we request a fiux density above 15 G on both footpoints. The resulting PDFs are shown 
in Fig. |5j The black and red lines correspond to the first and second time series, respec- 
tively. Additionally, we plot PDFs for distances between field line footpoints calculated by 
extrapolating calibrated data following Martinez Pillet et al. (2010), as explained in Sec. [2| 



which involves no inversion procedure (yellow and green curves, for the first and second time 
series, respectively). All the curves approximately fit an exponential distribution at scales 
between 1 and 10 Mm, i.e., at granular and mesogranular scales; the slope is rflog {PDF)/dx 
^ -0.6 Mm-^ Beyond 10 Mm, at supergranular scales, the curves show a deviation from 
the exponential. This bump most probably corresponds to the presence in that time series 
of strong network elements. Anyway, given the size of the IMaX field of view (~ 30 Mm) 
and the long duration of the largest convection cells, we cannot draw statistical inferences 
from the PDF for scales above, say, 10 Mm. On the other hand, we have carried out a 
Kolmogorov-Smirnov test of the goodness of fit of the footpoint distance data in the 1-10 
Mm range to a lognormal distribution, additionally to the exponential distribution; a log- 
normal distribution is found to fit well magnetogram data or, more generally, data resulting 



from the fragmentation of magnetic elements (Abramenko & Longcope 2005 see also Bog- 



dan et al. 1988). In our case, the test, carried out for individual snapshots in the given 
distance range, favors the exponential distribution. Details of this analysis will be given in 
a publication in preparation. 

We note the constancy of the slope of the PDF at scales between 1 and 10 Mm and that 
the characteristic decay distance is approximately 1.7 Mm. This is probably a consequence 
of the fact that there are no intrinsic horizontal scales in that range other than granulation, 
e.g., because mesogranulation is the direct result of other convection scales, rather than 
representing a primary scale in which energy is being injected into the convective fiows. 



6. Discussion and Conclusions 

An important open question in solar physics is the precise nature of the convection scales 
with size and duration above the granular values. While both granulation and supergranula- 
tion yield velocity field patterns that have been observed at the surface, mesogranules have 
been detected only through indirect proxies, like tracking of intensity patterns, which do not 
provide reliable evidence of underlying convection cells in that range of sizes and durations. 
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There is an ongoing debate on whether there is a continuum of sizes for the convection cells 
on scales above granular, possibly with self-similar properties and with no particular scale 



being singled out within that range (Nordlund et al. 2009), or whether the mesogranular 



scales are just the result of a collective interaction between families of granules. The mag- 
netic field can provide an alternative and more direct avenue to explore convective patterns 
since the magnetic flux can be measured directly using Stokes polarimetry techniques. The 
magnetic elements appear with a broad spectrum of flux densities ( Orozco Suarez et al.||2007 



Khomenko et al.|[2003 ) and can also be detected through G-Band or Ca-II bright points (see 



de Wijn et al. 2005 Sanchez Almeida et al. 2010) and therefore can provide important clues 



concerning the nature of the flows underlying the mesogranular scales as well as about the 
magnetic elements themselves. 

The high spatial and temporal cadences of the IMaX data allow us to try a few differ- 
ent, complementary studies of the relation between magnetic field and mesogranular flows. 
First, we have carried out Lagrange tracing of mass elements following horizontal flow fields 
obtained through LOT of intensity maps; we have obtained a number of improvements com- 
pared with traditional cork maps (faster development of the mesogranular lanes, no need for 
time averages). We have also obtained an upper bound for the width of the mesogranular 
lanes of some 280 km. Second, we have provided quantitative measures for the associa- 
tion between magnetic elements and mesogranular lanes. The large majority (85 %) of the 
magnetic elements with flux density above ~ 100 G are found in 120-km neighborhoods of 
mesogranular lanes with Pf,> 8 and about 80 % of flux elements above 30 G are located near 
mesogranular lanes with > 2. Our results indicate a good coupling between the flow field 
and the magnetic elements, suggesting that the evolution of the latter is mostly kinematic. 
Third, we have considered the connectivity between magnetic elements and, in particular, 
the distance between footpoints of field lines anchored in the photosphere. The probability 
density function of such distances shows that there is abundant connectivity on mesogranular 
scales; it also shows that the distribution is basically featureless in those scales, with only 
one characteristic value, the slope of the distribution, equal to (1.7 Mm)^^. 

Our results concerning statistics of separations of field line footpoints suggest that there 
is no intrinsic scale of convection in the mesogranular range. This may only mean that 
there is no mechanism for direct injection of energy into convection on those scales, but 
does not rule out the existence of convection cells with those sizes, e.g., through nonlinear 
interactions of cells at other scales (as in a turbulent cascade) or through the interaction of 



thermal downflows (Rast 2003). 
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Fig. 4. — Magnetic field lines showing the connectivity between flux concentrations. Back- 
ground image: map of (purple and blue: > 0; green and red: B^ < 0). The background 
10 X 10 Mm^ patch around the point (10,20) Mm in Fig. [l| 
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Fig. 5. — Probability density functions of the distance between footpoints of field lines. Black 
and yellow: data from the first time series (inverted data: black; data calibrated following 



Martinez Pillet et al. (2010): yellow). Red and green: same as black and yellow, but for 



series 2. The straight dashed line represents an exponential fit of the black curve for the first 
^10 Mm. 



